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1.  INTRODUCTION 

Any  attempt  to  go  past  the  limits  of  specific  energy  (energy  per  unit  mass)  imposed  by 
the  chemical  bond  in  the  production,  generation,  containment,  manipulation,  conversion 
or  utilization  of  energy  or  energetic  fluids  requires  the  application  of  the  principles  of 
magneto-  and  electrofluidmechanics.  As  examples,  we  may  consider  the  generation  of 
electricity  from  gases  at  temperatures  above  2000°  K,  the  acceleration  of  mass  to  velocities 
above  1600  m/sec  and  the  manipulation  of  materials  at  temperatures  above  3000°  K. 

Although  many  other  interesting  applications  of  magneto-  and  electrofluidmechanics 
exist  at  lower  temperatures  (e.g.,  pumping  of  liquid  metals)  all  these  other  applications 
are  technologically  less  demanding  by  comparison. 

There  is  not  much  doubt  that  this  is  a  very  broad  field  of  study  indeed,  since  it  involves 
all  of  classical  physics  and  thermodynamics  (through  the  coupling  of  the  Maxwell  equa¬ 
tions  with  the  laws  of  motion  and  conservation),  most  of  modern  physics  (e.g.,  through  the 
equation  of  state  and  atomic  and  molecular  ionization,  excitation,  and  radiative  processes) 
and  most  of  chemistry  (e.g.,  through  the  law  of  mass  action).  In  fact  few  fields,  if  any,  hold 

*  This  work  was  sponsored  in  part  by  the  U.S.  Air  Force  Office  of  Scientific 
Research  under  contracts  F49620-83-C-0115  anu  F49620-84-C-0068  under 
the  direction  of  Majors  Henry  Pugh  and  Bruce  Smith.  Their  unflinching 
support  under  adverse  conditions  is  gratefully  acknowledged. 
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more  promise  for  solving  some  of  the  outstanding  engineering  problems  of  mankind.  Great 
opportunities  exist  in  this  field  of  study  connected  with  such  fields  of  endeavor  as  edu¬ 
cation,  research,  and  development,  the  modernization  of  industry,  energy  self-sufficiency, 
and  others. 

Needless  to  say,  the  degree  of  complexity  of  this  field  of  study  is  commensurate  with 
the  number  of  disciplines  involved  and  in  fact  the  pace  of  success  in  it  is  intimately  coupled 
to  and  exposes  all  the  weaknesses  of  these  disciplines.  Without  careful  planning  there  is 
constant  danger  that  the  weakest  link  of  the  weakest  scientific  or  engineering  disciplines 
involved  will  dictate  the  pace  of  development  of  this  field. 

The  United  States  government  alone  has  spent  in  excess  of  8500  million  in  the  last 
ten  or  twenty  years  attempting  to  reduce  to  practice  or  commercialize  various  applications 
of  magneto-  and  electrofluid  mechanics.  The  magnitude  of  the  problems  and  the  original 
misunderstanding  exhibited  by  many  investigators  engaged  in  these  fields  were  such  that 
not  much  has  yet  materialized  from  all  these  efforts. 


2.  FUNDAMENTAL  THEORY 

2.1  The  Equations  of  Charged  and  Conducting  Fluids 

Magnetofluidmechanics  is  concerned  with  the  mechanics  of  fluids  which  may  support 
electromagnetic  body  forces.  These  forces  arise  from  the  presence  in  the  fluid  of  a  charge 
density  pc  and  a  current  density  J0.  The  total  current  density  J0  consists  of  a  convected 
portion  pc U  and  a  conduction  portion  J  which  is  carried  by  charge  drifting  at  a  relative 
velocity  with  respect  to  the  macroscopic  mean  velocity  U.  The  electromagnetic  force 
density  which  acts  upon  the  fluid  is  given  by 

f  =  pcE  +  J0  x  B  (2.1) 

where  E  and  B  are  the  electric  and  magnetic  fields  which  are  themselves  determined  by 
the  charge  and  current  present  both  in  the  fluid  and  its  surroundings.  The  corresponding 
power  flow  from  the  electromagnetic  field  to  the  fluid  is  given  by 

P  =  J0  •  E.  (2.2) 

The  macroscopic  (as  opposed  to  the  species)  conservation  laws  for  a  fluid  subject  to 
electromagnetic  forces  are 

~  +  pV  •  U  =  0,  (2.3a) 

p  =  ^ '  ^1  +  PcE  +  Jo  x  B,  (2.3 b) 
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p^(e  +  U2/2)  +  V  ■  (pU)  =  V  •  (n .  U)  +  V  •  q  +  Jo  •  E.  (2.3 c) 

The  mass  density,  velocity,  and  thermal  energy  per  unit  mass  of  the  fluid  are  p,U,  and  e. 
The  pressure  is  p,  II  is  the  viscous  stress  tensor,  and  q  the  heat  flux  vector. 

Since  each  species  (electrons,  ions  of  different  charge,  atoms  or  molecules  of  different 
mass)  in  the  fluid  may  have  its  own  temperature,  number  density,  etc.  and  its  own  con¬ 
stitutive  or  transport  properties,  Eqs.  (2.3)  are  more  properly  written  separately  for  each 
species  and  then  appropriately  summed.  However,  that  treatment  is  beyond  the  scope  of 
this  article.  It  can  be  found  in  such  sources  as  Demetriades  et  al.  1963,  1966,  1982  and 
Argyropoulos  et  al.,  1969,  and  Sutton  and  Sherman,  1964). 

In  most  engineering  applications  of  MHD  there  are  important  regions  where  the  elec¬ 
tron  temperature  is  not  the  same  as  the  gas  temperature  and  the  electron  number  density 
does  not  correspond  to  the  gas  temperature  through  some  appropriate  Saha  equation  (at 
the  gas  temperature)  but  is  closer  or  equal  to  the  number  density  obtained  from  a  Saha 
equation  at  the  electron  temperature  (Kerrebrock,  1962.)  it  is  important  to  retain  in  the 
system  of  Eqs.  (2.3)  at  least  the  separate  electron  energy  equation  which  couples  the  local 
electron  temperature  to  the  electric  and  magnetic  fields  and  local  electron  number  density. 

An  adequate  form  of  this  electron  energy  equation  is 

Je  •  E'  =  3/2 kneuSeff(Te  -  Ta)  +  R  +  ]T  i,c,  (2.3d) 

i 

where  Je  is  the  electron  current,  R  represents  the  radiative  losses  from  the  electrons,  is 
the  rate  of  “chemical”  reaction  of  electrons  and  heavy  particles  per  unit  volume  per  unit 
time,  e,  is  the  characteristic  energy  that  is  exchanged  in  each  such  reaction  i  in  units  of 
energy  per  reaction,  k  is  the  Boltzmann  constant,  ne  is  the  number  density  of  electrons,  v 
is  the  total  electron  collision  frequency,  Seff  is  the  effective  energy  loss  factor  per  electron 
collision  that  involves  elastic  and  inelastic  energy  exchange,  Tt  is  the  electron  temperature, 
Ta  is  the  heavy  species  (atoms,  molecules  or  ions)  temperature,  and  me  is  the  mass  of  the 
electron.  Usually,  for  diatomic  or  polyatomic  gases  and  at  temperatures  of  current  interest 
in  MHD,  R  and  s<  are  negligible  compared  to  the  collision  term  that  contains  ( Te  —  Ta). 

The  conservation  laws  Eqs.  (2.3)  must  be  supplemented  with  thermodynamic  state 
equations 


na  =  na(p,  T),  p  =  p(p,  T),  e  =  e(p ,  T).  (2.4) 

The  no,  are  the  number  densities  of  the  various  species  which  compose  the  plasma.  A 
preeminent  species  in  an  ionized  gas  is  the  electron  number  density  ne  whose  equlibrium 
state  equation  (2.4)  is  known  as  the  Saha  equation.  For  a  singly  ionizable  species  of  number 
density  n„  the  Saha  equation  has  the  simple  form 

r£  _  2nm,kTe  _(/kT, 
n.  ~  h2 
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where  h  is  Planck’s  constant  and  t  is  the  ionization  energy.  Although  we  have  indicated  no 
dependence  upon  the  electromagnetic  fields  in  the  state  relationships,  it  should  be  noted 
that  fluids  with  permanent  electric  or  magnetic  dipoles  possess  thermodynamic  property 
functions  which  depend  explicitly  upon  the  fields  E,  B  (de  Groot  and  Mazur,  1962;  Penfield 
and  Haus,  1967). 

The  constitutive  relationships  for  the  stress,  heat  flux,  and  conduction  current  consis¬ 
tent  with  the  Navier-Stokes  level  of  approximation  are 


(n)tJ  =  —2  pv 


dUA 

1W»-  1 

2  \dij 

dxi  J 

3dztSii, 

(2.5a) 


q  =  -AVT  +  t/>J, 


(2.56) 


J  =  a(E  +  U  x  B)  +  4>VT.  (2.5c) 

Equation  (2.5c)  for  the  current  J  which  complements  the  Navier  Stokes  relationship  for 
FI  and  q  in  a  conducting  fluid  is  called  the  “Ohm’s  Law.”  The  coefficients  p,  A,  a,  ip,  <p  are 
the  scalar  transport  coefficients  of  viscosity,  heat  conduction,  and  electrical  conductivity. 
In  general  the  heat  flow  and  conduction  current  are  coupled;  however  the  coupled  heat 
flux  term  <pVT  in  the  current  equation  may  be  neglected  whenever  the  electric  field  E'  = 
E  4-  U  X  B  is  much  greater  in  magnitude  than  (fc/e)VT  where  k/e  is  the  ratio  of  the 
Boltzmann  constant  to  the  electron  charge. 

When  a  magnetic  field  is  present,  the  transport  properties  become  tensorial.  This 
effect  is  known  as  the  Hall  effect  and  is  described  by  the  Hall  vector  0  =  (e/m)  B/i/  where 
e/m  is  the  charge  to  mass  ratio  of  the  conduction  carriers  and  u  is  their  effective  collision 
frequency  for  momentum  scattering.  When  0^1,  terms  proportional  to  VU  x  0,  VT  x  0, 
and  Ex/?  must  appear  in  Eqs.  (2.5a)  through  (2.5c)  respectively.  The  Hall  effect  on  the 
viscous  stresses  is  of  significance  in  fully  ionized  plasma  in  very  strong  magnetic  fields 
(Braginskii,195S;  Longmire  and  Rosenbluth,  1956).  In  partially  ionized  gases  the  Hall 
effect  upon  the  stress  tensor  is  generally  negligible  although  its  influence  upon  the  heat 
flux  and  current  may  be  substantial.  In  partially  ionized  gases  the  principal  contribution 
to  the  stress  tensor  arises  from  the  collisional  interaction  of  the  heavy  species  in  the  gas, 
most  of  which  are  not  charged. 

The  full  “generalized  Ohm’s  Law”  including  the  case  when  the  electron  temperature  is 
not  equal  to  the  heavy  species  temperatures  has  been  developed  (Demetriades,  Hamilton, 
Ziemer,  and  Lenn,  1962;  Demetriades  and  Argyropoulos,  1966;  Schweitzer  and  Mitchner, 
1967).  In  the  simplest  case  in  which  temperature  gradients  may  be  neglected,  Eq.  (2.5c) 
becomes 


Jl  = 


(1  +  02) 


(Ej.  -  E±  x  0 ), 


J||  =  <rE||, 


(2.7a) 

(2.76) 
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where  _L  and  ||  denote  directions  perpendicular  and  parallel  to  the  magnetic  field  respec¬ 
tively.  The  electric  field  E'  =  E  +  U  x  B  is  the  electric  field  in  a  frame  of  reference  moving 
with  the  fluid  velocity  U.  Using  Eqs.  (2.7)  the  expression  J  •  E'  may  be  expressed  as 

J-E'  =  72/"=(rf^)(£;1)2+^'- 

The  power  flow  J-E'  is  positive  definite  (“Joule  dissipation”)  and  represents  the  dissipation 
of  electrical  power  into  internal  thermal  energy  as  follows  from  the  thermal  energy  equation 
which  resuits  from  Eq.  (2.3c)  with  the  kinetic  energy  U2  / 2  eliminated  by  Eq.(2.3b). 

The  charge  density  pc  and  current  J0  determine  the  fields  E,  B  through  the  Maxwell 
equations: 


l/c23E fdt  —  V  x  B  =  —  poJo,  (2.8a) 

dB/dt  +  V  x  E  =  0,  (2.86) 

V-E  =  pc/co,  (2.8c) 

V  •  B  =  0.  (2.8d) 


Equations  (2.3)  and  (2.8)  with  the  state  and  constitutive  relationships  (2.4),  (2.5), 
and  (2.7)  form  a  closed  coupled  system  for  the  fluid  variables  p,  U,  e,  pc,  Jo  and  the  elec¬ 
tromagnetic  field  variables  E  and  B. 

2.2  Theorems  of  Bernoulli  and  Kelvin 

Using  the  identity  U  •  VU  =  VU2  /2  —  (V  x  U)  x  U,  the  momentum  equation  may 
be  expressed  as 


P(dv/dt -nxu  +  pVJ72/2)  +  v(P  +  w) 

=  -V  •  n  +  Po  lB  •  VB  -f  eo(E  •  VE  -  pcE),  (2'9) 

where  O  =  V  x  U  is  the  vorticity  and  w  =  (eo E2  4-  p~lB2)/ 2  is  the  electromagnetic  energy 
density.  When  the  right  hand  side  of  Eq.  (2.9)  vanishes  under  conditions  of  uniform  density, 
Bernoulli’s  theorem  for  electromagnetic  flow  is  obtained: 

U2/ 2  +  (e0E2  +  p;lB2)/2  +  p/p  =  C(t).  (2.10) 

Particular  velocity  fields  which  possess  Bernoulli  integrals  Eq.  (2.10)  include  all  veloc¬ 
ity  fields  which  may  be  represented  as  U  =  ki  x  E  -{-  k2  x  B  where  kj  and  k2  are  arbitrary 
vectors.  All  electromagnetic  fields  which  are  free  of  tangential  gradients  similarly  allow  the 
existence  of  a  Bernoulli  integral.  Most  general  flows  of  engineering  interest  do  not  satisfy 
these  conditions. 

The  curl  of  Eq.  (2.3b)  at  uniform  density  yields  the  vorticity  equation  in  electromag¬ 
netic  flow: 
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(2.11) 


P(Dn/Dt  -  n  •  vu)  =  -  v  x  (v  •  n) 

+  B  •  VJo  —  Jo  *  VB  +  peV  x  E. 

When  the  RHS  of  Eq.  (2.11)  vanishes  Kelvin’s  theorem  is  obtained:  the  flux  of  vorticity 
through  an  elemental  area  O  •  6  A  is  conserved.  In  conventional  flow,  viscous  stresses  act 
as  the  only  mechanism  for  the  creation  of  vorticity;  however  in  electromagnetic  flow,  the 
curl  of  the  Lorentz  force  does  not  generally  vanish  even  in  the  case  of  “ideaF  flow  in  which 
<r-1  — *  0  and  there  is  no  Joule  dissipation.  Gradients  of  the  current  density  in  the  direction 
of  the  magnetic  field  can  represent  an  important  source  of  vorticity  in  electromagnetic  flow. 

2.3  The  Magnetohydrodynamic  and  Electrohydrodynamic  Approximations 

The  divergence  of  Eq.  (2.8a)  and  use  of  Eq.  (2.8c)  with  Jo  =  pcU  +  J  yields  the 
equation  governing  the  charge  density: 


DpJDt  +  peV  •  U  =  —pc/e0.  (2.12) 

From  Eq.  (2.12)  rc,  the  decay  time  for  charge,  is  given  by  rc  =  e0/cr.  If  t  =  L/U  is 
a  characteristic  macroscopic  time  scale  where  L  is  the  length  scale,  the  condition  t/rc  = 
( L/tcU )  «  1  divides  conductors  from  non-conductors.  If  t/rc  1,  the  electrohydrodynamic 
approximation  J  «  0  holds  and 


DpJDt +  pcV-  U  =  0.  (2.13a) 

Correspondingly  the  Lorentz  force  becomes 

f  =  pc(E  +  Ux  B).  (2.136) 

For  t/rc  >  1  there  results  pc  «  0  and  Jo  =  J.  Equation  (2.12)  then  becomes 

V-J  =  0,  (2.14a) 

which  requires  that  Eq.  (2.8a)  takes  the  form 

V  x  B  =  ^0J.  (2.146) 

Correspondingly  the  Lorentz  force  becomes 

f  =  J  x  B.  (2.14c) 

The  approximation  embodied  in  Eqs.  (2.14)  which  stems  from  pc  =  0  is  known  as  the 
magnetohydrodynamic  approximation.  Note  particularly  that  the  neglect  of  dE/dt  from 
Eq.  (2.8b)  in  this  approximation  prohibits  a  description  of  electromagnetic  radiation. 
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2.4  The  Equation  of  Magnetic  Induction 

The  elimination  of  J  and  E  from  Eq.  (2.14b)  and  the  Ohm’s  Law  with  Eq.  (2.8b)  re¬ 
sults  in  the  induction  equation  for  B  for  uniform  density  and  electrical  transport  properties 
<r,  0  : 


DB/Dt  -  B  ■  VU  =  77V2B.  (2.16) 

Equation  (2.16)  describes  the  convection  and  diffusion  of  the  magnetic  field  with  the 
magnetic  diffusivity  rj  =  l//x0cr.  When  r?  — ►  0.  the  electrical  forces  are  non-dissipative  and 
the  magnetic  field  is  convected  as  if  the  lines  of  force  were  bound  to  the  fluid.  The  left 
hand  side  of  Eq.  (2.16)  is  formally  identical  to  the  vorticity  equation  Eq.  (2.11):  B  •  6 A. 
the  flux  of  B  through  an  elemental  area  6 A,  is  conserved  throughout  the  flow  just  as  the 
flux  of  vorticity  O  ■  8 A  is  conserved  when  Kelvin’s  theorem  holds. 

The  time  scale  rg  for  the  field  to  diffuse  over  a  length  scale  L  is  found  from  (2.16)  to 
be 


rB  =  L2/ho<t-  (2.17) 

The  convective  time  scale  L/U  may  be  compared  to  the  diffusion  time  scale  rg  in  a  non- 
dimensional  parameter  known  as  the  magnetic  Reynolds  number  Rm  (Lundquist,  1952): 


Rm  =  tb/(L/U)  =  WUL.  (2.18) 

For  Rm  -C  1,  magnetic  fields  diffuse  rapidly  through  the  region  of  interest.  The  regime 
of  vanishing  Rm  has  other  important  consequences.  Consider  a  steady  magnetic  field  Bo 
which  is  established  by  currents  external  to  the  fluid.  Let  J  be  the  current  present  in  the 
fluid  and  B'  given  by  Eq.  (2.14b)  the  field  induced  by  the  currents  in  the  fluid.  From 
Ohm’s  Law  with  E  «  U  x  (B0  +  B')  the  current  is  given  by  J  «  cr\J  x  (B0  +  B').  The 
field  B'  induced  by  this  current  is  from  Eq.  (2.14b) 


VxB'«  no  J  «  A*o^U  x  (B0  +  B').  (2.19) 

Since  |V  x  B'|  ss  B'/L,  there  results 

5^-57  «  wUL  =  R„.  (2.20) 

When  Rm  >  1,  the  field  satisfies  B'/Bo  «  1;  but  for  Rm  <.  1,  this  condition  becomes 
B'/Bq  «  Rm  1  and  the  induced  field  B'  may  be  neglected  compared  to  the  applied  field 
B0.  From  Eq.  (2.8b)  we  have  dB/dt  w  —V  x  E.  For  E  ss  U  x  B'  and  dfdt  of  the  order 
of  the  convective  time  scale  U fL  it  will  be  true  that 


\dB' fdt\  ~ 

TvTeT ' 


(2.21) 
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In  the  approximation  of  low  magnetic  Reynolds  number  Ohm’s  Law  is  expressed  only  in 
terms  of  applied  field  Bq: 


J=  — - 1— -|E  +  UxBo-Ax(E  +  UxBo)],  (2.22) 

( 1  +  a0 / 

while  the  Maxwell  equation  (2.8b)  becomes 


V  x  E  =  0,  (2.23) 

and  Eq.  (2.8a)  in  the  MHD  approximation  becomes 

VJ  =  0.  (2.24) 

Equations  (2.22)-(2.24)  serve  to  completely  determine  the  fields  E,  J  in  the  low  mag¬ 
netic  Reynolds  number  and  MHD  approximations.  The  magnetic  field  B  C  Bo  induced  by 
the  current  governed  by  Eqs.  (2.22)  -  (2.24)  may  be  determined  from  Eq.  (2.14b). 

2.5  Characteristic  Parameters 

The  characteristic  parameters  for  the  general  fluid  equations  in  absence  of  body  forces 
are  three:  The  Mach  number  M  =  U/c  (where  c  =  yjjp/ p  is  the  speed  of  sound),  the 
viscous  Reynolds  number  Re  =  ULfv ,  and  the  Prandtl  number  Pr  =  Cpv/ p\  where  Cp  is 
the  specific  heat  at  constant  pressure.  Additional  parameters  are  required  in  magnetohy¬ 
drodynamic  flow. 

The  first  of  these  parameters  is  the  magnetic  Reynolds  number  Rm  which  describes 
the  convection  of  magnetic  energy  compared  to  its  diffusion.  One  other  independent  one¬ 
dimensional  number  is  required  for  the  description  of  electromagnetic  effects  in  the  MHD 
approximation.  For  Rm  3>  1  the  ratio  of  magnetic  energy  density  to  the  fluid  energy 
density  is  an  appropriate  choice: 


lu  =  B2/p0pU2.  (2.25) 

A  second  choice  is  the  ratio  of  magnetic  pressure  to  fluid  pressure: 

Ip  =  B2/p0p.  (2.26) 

The  parameters  Iu,  Ip  are  not  independent  but  are  related  through  the  Mach  number  M  : 

Ip  =  7  M2Iu-  (2-27) 

Under  conditions  where  Rm  C  1,  magnetic  diffusion  effects  dominate  over  convection 
and  a  more  appropriate  choice  of  parameters  are  the  interaction  parameters  Su  and  Sp. 
The  parameter  Su  is  defined  as  the  ratio  of  the  Lorentz  force  (with  J  expressed  from  Ohm’s 
Law)  to  the  change  of  fluid  momentum  in  length  L : 

Su  =  <rB2L/pU.  (2.28) 
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The  parameter  Sp  is  the  ratio  of  the  Lorentz  force  (similarly  expressed  through  the  Ohm’s 
Law)  to  the  pressure  gradient  over  length  L  : 

Sp  =  aUB2Lfp.  (2.29) 

The  parameters  Su-,  Sp  are  not  independent  of  Iu,Ip  hut  are  related  through  the  magnetic 
Reynolds  number: 


Su  =  Rmlu,  (2.30a) 

Sp  =  RmIv.  (2.306) 

The  Hartmann  number  Ha  is  defined  as 

Ha  =  y/'SuRe-  (2.31) 

and  is  a  measure  of  the  ratio  of  the  Lorentz  force  to  the  viscous  force. 

For  general  magnetohydrodynamic  flow,  a  typical  and  complete  set  of  independent 
parameters  are  il /,  Re,  Pr ,  i?m,  and  Su-  When  the  electrical  conductivity  is  tensorial  due 
to  Hall  effect,  the  Hall  parameter  8  must  be  added  to  this  set. 


3.  BASIC  MAGNETOHYDRODYNAMIC  PROCESSES 


3.1  Local  Currents  and  Fields 


Consider  the  nature  of  the  currents  and  fields  J,E  in  a  local  region  of  flow.  From  the 
Ohm’s  Law  Eq.  (2.7)  it  can  be  seen  that  the  magnetic  field  effects  involve  U  x  B  and 
(U  x  B)  x  8.  These  effects  are  confined  to  the  plane  perpendicular  to  B.  Orient  B  in  the  z 
direction  so  that  B  =  B(0,0,  B),  U  in  the  x  (axial)  direction  with  U  =  U(17,0,0),  and  let 
E  =  E (Ex,  Ey,0)  where  y  represents  the  transverse  direction.  The  axial  and  transverse  load 
factors  Kx,Ky  which  serve  to  summarize  the  nature  of  the  exterr  al  electrical  connection 
to  the  conducting  fluid  are  defined  as 

K  -Et 

T  3UB{\-Ky)' 

(3.1a) 

Ky  =  Ey/UB. 

(3.16) 

The  Ohm’s  Law  Eqs.  (2.7)  then  becomes 

Jz  =  (l-Kx)(l-Kv)^^r)8UB, 

(3.2a) 

J^-il-KyKl  +  K^jf^UB. 

(3.26) 
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For  Ky  =  1  both  the  axial  and  transverse  currents  are  “open  circuited”  and  from  Eqs. 

(3.1)  the  open  circuit  values  of  the  fields  are  Ez  =  0  and  Ev  =  UB.  For  Kz  =  1  for  any 
Ky,  the  axial  current  vanishes  and  the  transverse  current  is  given  by 

Jy  =  -(1  ~Kv)aUB.  (3.3) 

Thus,  if  Jz  is  prevented  from  flowing  the  transverse  current  Jy  achieves  a  ievel  equal  to 
that  which  it  would  have  in  the  absence  of  the  Hall  effect;  an  axial  or  ’’Hall”  field  Ez  will 
be  present,  however,  given  by  Ez  —  —BEy  where  Ey  =  (1  —  Ky)UB.  Clearly,  for  large 
Hall  parameter  tf  the  axial  field  will  dominate  the  transverse  field.  On  the  other  hand  if 
A'r  =  0,  the  axial  field  is  shorted  and  the  transverse  current  from  Eqs.  (3.2)  is  reduced  by 
the  factor  1/(1  +  02)  from  its  value  in  the  unshorted  case  given  by  Eq.  (3.3). 

3.2  Force,  Power,  and  Dissipation 

The  local  Lorentz  force  f(/*,/v,0)  and  power  P  are  given  from  Eqs.  (2.1),  (2.2)  and 

(3.2)  as 


/*  =  -(1  —  Ky)(l  +  KXP2) 


(1  +  B2) 


UB 2, 


fv  =  “(I  ~Kx){l~Ky) 


d  +  02) 


BUB2, 


<jU2B2 


P  =  -(1  -  Ky)  [AVI  +  KrB2)  +  (1  -  A'y)(l  -  Kx)Kxp2} 


(3.4a) 

(3.46) 

(3.5a) 


The  dissipated  power  P'  =  J  •  E'  is  given  by 


P'  =  (l- KyY(l+02Kl) 


aU2B2 


(3.56) 


For  Kz  <  1,  Kv  <  1  the  force  fx  is  opposed  to  the  flow  velocity  U  and  power  (P  <  0) 
is  extracted  from  the  flow.  For  Ky  >  1.  the  force  fx  is  in  the  direction  of  the  flow  and 
power(P  >  0)  must  be  supplied  to  the  fluid  by  the  fields.  For  Kx  >  0,  a  transverse  force 
fy  will  exist  in  addition  to  the  axial  force  fx  if  0  >  0.  For  strong  Hall  effect  i?  >  1  this 
force  will  dominate  the  axial  force  fx  creating  significant  transverse  pressure  gradients  and, 
under  some  circumstances,  secondary  flow. 

The  total  axial  force  Fz  and  power  P  which  can  be  developed  in  a  volume  V  are  given 

by 


Fx  =  -  [  dv{  1  -  Ky)aUB2, 

Jv 

V=~  l  dv(\  -  Ky)KvcrU2B2. 
Jv 


(3.6) 

(3.7) 


Handbook  of  Fluid  Dynamics 


-  10  - 


©1987  STD  Research  Corporation 


In  a  uniform  system  these  expressions  have  the  simple  forms 


Fx  =  -(1  -  Ky)oUB2V,  (3.8) 

V  =  —(1  —  Ky)KyaU2B2V.  (3.9) 

It  is  important  to  observe  that  Eqs.  (3.8)  and  (3.9)  follow  from  Eqs.  (3.6)  and  (3.7) 
only  in  the  case  that  Ky,<7,l\  and  B  are  uniform  over  the  volume  V.  In  most  realistic 
magnetohydrodynamic  flows  these  quantities  are  significantly  non-uniform  (Demetriades, 
Argyropoulos,  and  Kendig,  1967:  Demetriades  and  Argyropoulos,  1969)  and  the  total  force 
or  the  total  power  cannot  be  reliably  calculated  from  Eqs.  (3.8)  and  (3.9). 


3.3  Magnetic  Meter  and  Pump 

Consider  a  channel  aligned  in  the  x  direction  with  crossed  electric  E(0,Ev,0)  and 
magnetic  B(0,0,i?)  fields.  The  steady  one-dimensional  incompressible  flow  equations  in 
the  low  magnetic  Reynolds  number  limit  with  /*  given  by  Eq.  (3.4a)  are 


(3.10) 


-^=  -(1  -Ky)oUB2.  (3.11) 

where  p0  =  p  +  pU 2  /2  is  the  stagnation  or  total  pressure.  Since  the  condition  V  x  E  =  0 
holds  we  have  Ey  =  -$/h  where  $  is  the  voltage  across  the  gap  h.  From  Eq.  (3.1  b)  this 
voltage  is  given  by  $  =  —KyUBh.  The  pressure  change  Apo  in  length  L  of  the  channel  is 
found  from  Eq.  (3.11): 


Apo  =  -(1  -  Ky)aUB2L.  (3.12) 

The  total  current  per  unit  depth  w  which  flows  in  length  L  is  given  by 

I/w  =  JyL  =  -(1  -  Ky)aUB2L.  (3.13) 

The  pressure  change  can  thus  be  represented  as 

Apo  =  (I/w)B.  (3.14) 

The  power  delivered  to  the  fluid  per  unit  depth  is 

PL  =  JyEyL  =  -(I/w)<S>.  (3.15) 

At  open  circuit  1  =  0  the  voltage  at  this  condition  is  given  by  $  =  $oc  =  —UBh.  This 
voltage  can,  in  principle,  be  used  to  meter  the  average  velocity  U  which  exists  in  the  channel 
when  a  known  magnetic  field  is  applied  over  the  gap  of  known  dimension  h.  In  practice 
such  meters  must  be  corrected  for  leakage  current  losses  and  polarization  voltages  which 
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exist  on  the  electrodes  for  small  devices.  For  a  device  of  dimension  h  =  10cm  and  imposed 
magnetic  field  B  =  0.1  Tesla,  a  velocity  of  10 m/s  will  register  only  lOOmV.  Faraday 
attempted  to  measure  the  velocity  of  the  Thames  River  using  this  principle.  Although 
h  «  100m  for  the  river,  the  earth’s  magnetic  field  is  only  of  the  order  \0~*Tesla  and 
thus  $oc  «  10m V'  for  U  «  1  m/s  which  is  in  the  noise  range  of  the  spurious  leakage  and 
polarization  effects. 

To  operate  the  channel  as  a  pump  (Apo  >  0)  a  voltage  $  <  —UBh  must  be  imposed 
across  the  gap.  At  the  no-flow  condition  U  =  0,  the  pressure  rise  is  maximum  and  is  given 
by 


(Apo  )max  =  (-$)cBL/h.  (3.16) 

For  liquid  metals,  the  pumping  voltage  is  low  because  of  the  high  electrical  conductivity. 
For  L/h  «  1,  B  «  1  Tesla,  a  10 7 S/m,  and  (Apo)max  »  l atm.  the  pumping  voltage  4> 
is  of  the  order  $  «  1  volt.  The  corresponding  current  from  Eq.  (3.13)  is  (//tu)  «  10 7 A/m 
and  the  power  per  unit  depth  is  Pi  =  100 kW/m.  These  voltage-current  characteristics 
place  a  severe  demand  on  the  pump  power  supply. 


4.  FUNDAMENTAL  MAGNETOHYDRODYNAMIC  FLOWS 


4.1  Quasi  One-Dimensional  Compressible  Flow 

In  the  quasi-one-dimensional  description,  the  flow  variables  are  averaged  quantities 
over  the  cross-  section  of  a  duct.  The  conservation  laws  Eqs.  (2.2)  integrated  over  the 
cross-section  with  the  variables  now  representing  averages  over  the  cross-section  A  are 


d(pA)  d(pUA) 
dt  +  dx 


d(pUA)  d(pU2A) 
dt  dx 


~  fx)  -  TwC  4-  Ey, 
ox 


(4.1  a) 
(4.16) 


d[p(U2/2  +  e)A]  d(pUhpA) 
dt  dx 


PA  4-  4"  Et. 


(4.1c) 


The  dependent  fluid  variables  are  the  pressure  p,  enthalpy  h  =  e  4-  p/p,  and  stagnation 
enthalpy  h0  =  h  +  U2  /2.  The  local  duct  cross-sectional  area  is  A  and  the  local  perimeter  is 
C.  The  functions  Ey,  and  E«  describe  the  sources  of  mass,  momentum  and  energy  due 
to  mass  addition  (or  depletion)  within  the  flow  train.  In  addition  to  the  conservation  laws, 
the  caloric  and  kinetic  equations  of  state  (2.3)  complete  the  description.  The  average  wall 
shear  stress  over  the  cross-section  is  tw  and  the  average  heat  flux  is  qw.  These  averages 
are  composites  of  the  electrode  wall  and  insulating  wall  shear  stresses  and  heat  fluxes. 
These  quantities  may  be  computed  from  the  boundary  layer  equations  for  each  wall.  The 
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results  of  these  boundary  layer  calculations  may  be  summarized  in  the  local,  dynamic 
friction  factor  C/(x,  c,  t)  and  Stanton  number  St(x,  c,  t)  where  c  is  a  perimeter  coordinate. 
The  wall  shear  and  heat  flux  are  directly  computed  in  terms  of  these  quantities  and  then 
averaged  over  the  cross  section  to  obtain  the  average  wail  stress  and  heat  flux  which  appear 
in  Eqs.  (4.1)  (Oliver.  Crouse,  Maxwell  and  Demetriades.  1980). 

The  most  fundamental  solutions  of  Eqs.  (4.1)  which  illustrate  electromagnetic  effects 
are  those  for  steady,  uniform  area  flow  in  the  absence  of  viscous  and  heat  conduction 
effects.  Eqs.  (4.1)  are  three  simultaneous  equations  for  the  primitive  variables  p,U,e. We 
may  eliminate  these  variables  in  favor  of  one  equation  in  the  Mach  number  M  =  U /  \J 7 p/  p 
using  Eqs.  (3.4a)  and  (3.5a)  for  fx  and  P  with  Kx  =  1  : 


1  dM 2 
\P  dx 


(7— l)(l+7  M2) 

v  7(1  -m 


[(l-A-y)(A\-A-v)], 


where  A', (A/)  is  defined  as 


(4.2) 


A'.(M) 


__v _ 1  +  (7  -  1  )/2  M2 

( 7  —  1 )  1  +  7  A/2 


In  Eq.  (4.2),  SJ,  =  aUB2 /p  is  the  pressure  interaction  parameter  per  unit  length.  This 
equation  defines  the  conditions  which  drive  the  flow  to  or  from  the  choke  point  M  —  1.  For 
all  Ky  <  1  the  Lorentz  force  opposes  the  flow  and  drives  the  flow  to  the  sonic  point.  For 
all  Ky  >  1  the  Lorentz  force  acts  in  the  direction  of  the  flow  and  drives  the  flow  away  from 
the  sonic  point.  When  (1  —  Ky)  %  1  holds  the  effect  of  the  Joule  heating  term  is  always  to 
drive  the  flow  to  the  sonic  point:  there  are,  however,  values  of  Ky  for  which  the  bracketed 
term  on  the  right  hand  side  of  Eq.  (4.2)  vanishes  and  continuous  passage  through  M  =  1 
in  an  “electromagnetic  throat”  is,  in  principle,  possible  (Neuringer.  1963).  It  can  be  seen 
that  there  are  two  possible  “throats”  at  Ky  =  1  and  Ky  —  A*(l)  =  7/2(7  —  1). 

Exact  integrals  for  steady  flow  in  the  absence  of  friction  and  heat  transfer  exist  (Rosa, 
1968,  pp.  51-54).  For  the  general  case  of  flow  which  is  unsteady  and  in  which  viscous  and 
heat  conduction  effects  are  important  the  solutions  are  best  obtained  by  finite  difference 
methods.  Such  solutions  of  Eqs.  (4.1)  are  presented  in  Section  6. 


4.2  Viscous,  Incompressible,  Magnetohydrodynamic  Flow 

Equations  (2.3)  and  (2.5a)  for  steady,  incompressible,  low  Rm  flow  in  the  x  direction 
between  the  planes  z  =  ±h  in  the  presence  of  B  =  B(0,0,I?),  U  =  U(l/x,0, 0),  J  = 
J(0,  Jy,0)  and  E  =  E(0,  Ey,Q)  reduce  to 


(4.3a) 


dp  r  o  d2Ut 

-=j,b+pv— 


(4.36) 
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Jy  =  —&(Ey  —  UXB).  (4.3c) 

Since  Ux  =  Ur{z),  define  Ky  as  Ky  =  Ey/(UX)B  where  (Ux)  is  the  average  velocity  between 
the  planes  s  =  ±h.  The  Ohm’s  Law  Eq.  (4.3c)  then  becomes 


Jy  =  -a(Uz-(Ux)Ky)B.  (4.4) 

It  can  be  seen  that  the  current  and  the  Lorentz  force  reverse  at  the  point  where 
Uz  =  Ky(Uz)  with  Jy  >  0  for  Ux  <  KV(UZ)  and  Jy  <  0  for  Uz  >  Ky(Uz).  Combining 
Eqs.  (3.3b)  and  (3.4)  and  observing  that  d2p/dxdz  =  0  since  there  is  no  motion  in  the 
transverse  direction  there  results 

■j/r  --j$-[Vz-  (K,  +  s;'){v,)]  =  o,  (4.5) 

where  Ha  =  Bh\Jcr / pv  is  the  Hartmann  number  and  Sv  =  aB2{Uz)/(—dp/dx) is  the 
pressure  interaction  parameter. 

The  solution  to  Eq.  (4.5)  for  Uz(z)  satisfying  the  boundary  conditions  Uz{h)  = 
Uz(—h)  is 

v,(z)  •  (K,  +  s;')(v,) 

By  taking  the  average  of  Eq.  (4.6),  the  interaction  parameter  Sp  may  be  represented  as  a 
function  of  Hartmann  number  and  load  factor: 


1  - 


cosh  Hazjh 


cosh  H„ 


(4.6) 


a&jUJ  g.-taph  H, 

'  l-dp/dx)  (1  +  A.  tan h  Ha  11 

Equation  (4.7)  may  be  used  to  determine  the  average  velocity  in  terms  of  the  Hart¬ 
mann  number,  load  factor,  and  the  pressure  gradient.  For  Ha^>  1,  we  find  Sp  «  1/(1  — fif,) 
with  the  Lorentz  force  balancing  the  pressure  gradient.  For  Ha  <£  1 ,  we  have  tanh  Ha  « 
Ha  —  Ha/3  and  Sp  «  H2  so  that  the  pressure  gradient  is  balanced  against  the  viscous 
force.  Using  Eq.  (4.7),  the  velocity  profile  (4.6)  may  be  expressed  edxclusively  in  terms  of 
Hartmann  number: 


Ur/{UX) 


Ha 

Ha  —  tanh  Ha 


1 


cosh  Haz/h 
cosh  Ha 


(4.8) 


The  velocity  profiles  for  insulating  walls  are  shown  in  Fig.  4.1  as  a  function  of  Hartmann 
number.  The  effect  of  the  magnetic  field  is  to  flatten  the  velocity  profile  in  the  core  where 
Jv  <  0,  /*  <  0  and  to  enhance  the  profile  near  the  sidewalls  where  Jy>  0  and  fx  >  0. 


5.  LARGE  MAGNETIC  REYNOLDS  NUMBER 
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Flow  at  large  magnetic  Reynolds  number  is  characterized  by  the  presence  of  magnetic 
fields  induced  by  currents  flowing  in  the  fluid  which  are  comparable  to  externally  imposed 
magnetic  fields.  These  internal  fields,  set  up  by  the  internal  currents,  tend  to  oppose  the 
imposed  field  (Lenz’  Law)  thereby  driving  the  total  magnetic  field  from  the  conducting 
fluid.  At  large  magnetic  Reynolds  numbers  the  currents  and  fields  are  confined  to  a  thin 
zone  at  the  boundaries  of  the  conducting  fluid. 


5.1  Magnetic  Field  of  a  Conducting  Fluid  Slab 

A  fluid  with  a  uniform  velocity  U  =  U(t/j,0,0)  moves  in  a  channel  of  dimension 
a  <C  L  where  Lis  the  scale  of  variation  of  an  imposed  magnetic  field  B0  =  Bo(0,0, B0)- 
Scalar  conductivity  of  magnitude  <r  is  switched  on  between  the  planes  x  =  O.a.  The 
conductivity  vanishes  outside  this  slab. 

The  governing  equation  for  Bz(x)  in  steady  flow  is  Eq.  (2.16)  (with  0  —  0): 


„SB,  92B, 
1  dx  '  dx*  ’ 

while  the  current  is  determined  from  Eq.  (2.14b): 

dBz 


jy  =  HO ' 


dx 


(5.1) 


(5.2) 


The  boundary  values  are  Bz(a),Bz{0).  The  total  current  per  unit  depth  flowing  through 
the  slab  is  given  in  terms  of  the  boundary  values  by 


I 


(5.3) 


This  current  is  also  given  in  terms  of  the  load  factor  Ky  =  Ey/UBo  from  Eq.  (3.3)  as 


=  -/' 


crUBo{Bz/Bo  -  Ky)dx. 


(5.4) 


The  boundary  conditions  to  be  applied  to  Eq.  (5.1)  is  that  the  induced  magnetic 
field  immediately  in  front  of  the  plasma  must  exactly  oppose  the  corresponding  induced 
magnetic  field  at  the  trailing  edge  of  the  plasma  so  as  to  maintain  the  condition  V  •  B  =  0. 
This  is  equivalent  to  demanding  that  the  average  of  the  boundary  values  be  equal  to  the 
imposed  field  Bq:  Bz(a)  +  B*(0)  =  2Bo- 

The  solutions  to  Eqs.  (5.1)  and  (5.2)  subject  to  these  conditions  (Demetriades  st.  al., 
1985)  are 


and 


Bz  =  BoKy  + 


2B0(1  —  Ky)  B„x/tt 

(1  +  e*m) 


(5.5) 
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(5.6) 


2BooUt(\  -  K„)eR'*x'a 
(1  +  eRm ) 


where  Rm  =  noUaa  is  the  magnetic  Reynolds  number.  From  Eqs.  (5.5)  and  (5.6)  it  can 
be  seen  that  the  field  and  current  are  progressively  confined  to  a  thin  zone  at  the  front  of 
the  conducting  slab  as  the  magnetic  Reynolds  number  is  increased  (Fig.  5.1).  The  total 
power  extracted  from  the  conducting  slab  per  unit  of  cross-section  flow  area  A  may  be 
expressed  as 


JyEydx  =  -4A'y(l  —  A'y ) ( B\ /2 ) fio ) U  tanh(Am/2). 


(5.7) 


As  the  magnetic  Reynolds  number  is  increased  we  have  tanh( Rm/2)  —*  1  and  the  extracted 
power  is  limited  to  the  product  of  the  magnetic  energy  of  the  applied  field  Bo  and  the  flow 
velocity  U. 


5.2  Planar  Flow  in  an  Inhomogeneous  Magnetic  Field 

The  behavior  of  the  fields  and  currents  in  flow  through  an  imposed  localized  magnetic 
field  at  large  magnetic  Reynolds  number  is  inherently  multi-  dimensional.  The  imposed 
field  Bo  =  Bo(0, 0,  Bo)  may  be  represented  as 


B0(x)  =  B./(x),  (5.7) 

where  Bm  is  a  constant  amplitude  and  f{x)  is  a  shape  function  with  a  maximum  at  /(0). 
The  fluid  is  contained  within  the  planes  y  =  ±h  and  flows  steadily  at  a  uniform  velocity 
U  =  U(t/r,0,0).  Conductors  are  arranged  on  the  planes  y  =  ±h  over  the  region  jx|  <  a. 
The  region  |x|  >  a  is  an  insulating  surface. 

The  governing  equation  for  Bz(x,y)  is  the  two-  dimensional  form  of  Eqs.  (2.16): 


d2Bz\ 
dy 2  )' 


and  the  current  J  =  J(JX,  7^,0)  is  given  by 


(5.8) 


dBx 
dy  ’ 


(5.9) 


The  boundary  conditions  are  that  the  current  vanish  at  x  =  ±oo  and  hence  dBz/dx  =  0 
at  x  =  ±oo.  On  the  conductors  the  tangential  field  Et  =  Jxl<r  must  vanish  while  the 
normal  current  must  vanish  on  the  insulating  surfaces,  i.e.,  dBz/dy  =  0. 

Because  of  the  nonuniform  shape  /(x)  and  the  mixed  boundary  conditions  on  y  =  ±h 
the  solutions  of  Eqs.  (5.7)-(5.10)  are  obtained  using  finite-difference  methods  (Oliver, 
Swean,  and  Markham,  1981).  The  appropriate  non-dimensional  parameters  are  the  mag¬ 
netic  Reynolds  number  Rm  =  \iooU a  and  the  current  load  parameter  ft0I/wB,  where  the 
total  current  per  unit  depth  flowing  through  the  conductors  to  the  external  circuit  is  I/w. 


Handbook  of  Fluid  Dynamics 


-  16  - 


©1987  STD  Research  Corporation 


In  Fig.  5.2  the  open  circuit  and  loaded  condition  for  a  magnetic  Reynolds  number 
Rm  =  1  are  exhibited.  It  can  be  seen  that  there  is  an  onset  of  eddy  current  cells  induced 
in  the  regions  of  gradients  in  the  imposed  magnetic  field  and  a  convection  of  the  current 
pattern  downstream  of  the  conductors.  The  principal  current  flowing  in  the  conductors, 
however,  is  still  confined  to  the  region  between  ti  <•  electrodes. 

In  Fig.  5.3  the  same  situation  is  represented  but  for  a  magnetic  Reynolds  number 
Rm  =  7.  Considerable  current  flow  induced  by  the  magnetic  field  gradients  exists  under 
open  circuit  conditions.  Note  how  the  upstream  eddy  cell  couples  into  the  electrodes. 
Under  load,  the  current  exists  completely  downstream  from  the  electrodes 


6.  MAGNETOHYDRODYNAMIC  POWER  GENERATORS 

Magnetohydrodynamic  power  generators  are  compact,  high  energy  density  converters 
of  thermal  energy  to  electricity.  Such  generators  may  function  as  “electro  magnetic  tur¬ 
bines”  in  a  complete  power  system  cycle  with  either  a  fossil  or  nuclear  heat  source  and  a 
steam  bottoming  plant.  In  this  function  they  serve  as  high  temperature  turbines  capable 
of  handling  a  top  cycle  temperature  in  excess  of  2500° K  as  weil  as  siag-laden  working  fluid 
in  the  case  of  coal  fired  plants.  Overall  MHD  topped  plant  efficiencies  should  be  in  the 
range  of  50  —  60%  with  the  MHD  generator  extracting  20  —  25%  of  the  input  chemical 
energy  (Heywood  and  Womack,  1969).  MHD  generators  may  also  function  as  portable 
devices  with  a  rocket  type  combustor  for  short  bursts  of  power  (Maxwell  and  Demetriades. 
1986;  Velikhov.  1975);  or  they  may  be  inherently  unsteady  driven  by  high  energy  explosives 
(Bangerter,  Hopkins,  and  Brogan,  1975). 

6.1  Ideal  Generators 

The  fundamental  starting  point  of  low  magnetic  Reynolds  number  generator  designs 
are  the  expressions  for  force  and  power  in  a  basic  crossed-field  situation,  Eqs.  (3.4)  and 
(3.5).  Three  principal  generator  types  follow  from  these  expressions:  Faraday,  Hall,  and 
Diagonal.  In  the  Faraday  generator  (Fig.  6.1a)  the  electrode  walls  must  be  segmented 
and  multiple  loads  attached  to  preserve  the  condition  of  vanishing  Hall  current  ( Jx  =  0, 
Kx  =  1).  The  load  current  in  axial  length  L  is  given  by  I  =  J9Lw  and  the  voltage  across 
the  electrodes  of  separation  h  by  $  =  —Evh.  This  configuration  possesses  only  an  axial 
Lorentz  force  and  has  a  power  density  which  is  not  degraded  by  the  Hall  effect.  The  large 
number  of  load  circuits  are  a  disadvantage  in  this  configuration.  In  the  Hall  generator 
(Fig,  5.1b)  the  Faraday  circuits  are  shorted  ( I\9  =  0)  and  tho  ends  of  the  machine  are 
connected  to  a  single  load.  The  load  current  for  a  machine  of  length  L  and  transverse  area 
hw  becomes  I  =  Jthw  and  the  load  voltage  for  this  configuration  is  $  =  -EZL.  This 
machine  is  only  efficient  at  large  Hall  parameter  /?;  in  addition,  a  large  transverse  Lorentz 
force  f9  exists  which  exceeds  the  axial  Lorentz  force  and  develops  secondary  flows  and 
transverse  pressure  gradients  which  degrade  performance. 
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The  diagonal  generator  is  constructed  much  like  a  Faraday  generator.  From  Eqs.  (3.1 ) 
the  angle  a  which  the  equipotential  lines  make  which  are  orthogonal  to  field  lines  is  given 

by 


tan  a  = 


2(1 -A'  )AV 


(6.1) 


For  a  given  Hall  parameter  and  operating  load  condition  Ky,Kx  the  equipoten- 
tials  which  are  orthogonal  to  the  field  lines  given  by  Eq.  (6.1)  may  be  made  of  con¬ 
ductors  and  the  power  again  extracted  by  a  single  load  attached  to  the  diagonal  elec¬ 
trodes.  This  machine  has  the  single  load  advantage  of  the  Hall  machine  without  its  per¬ 
formance  degradation.  In  this  case  the  load  current  into  a  diagonal  conductor  is  given  by 
I  =  hw(Jt  —  tan aJy)  and  the  load  voltage  by  $  =  —EZL. 

The  local  electrical  efficiency  T)  is  defined  as  the  ratio  of  the  electrical  power  P  to  the 
mechanical  power  flow  f  •  U.  The  force,  power,  and  efficiency  of  these  configurations  are 
given  in  Fig.  6.2.  It  will  be  noted  that  the  diagonal  configuration  becomes  inefficient  at 
operating  conditions  away  from  its  design  point. 

The  maximum  output  power  density  for  these  machines  is  aU2B2/ 4.  For  fossil  fuel 
combustion  products  seeded  with  an  easily  ionized  substance  such  as  a  potassium  or  cesium 
salt,  the  achievable  conductivity  is  of  the  order  of  10 S/m,  the  velocity  of  order  103m/s, 
and  the  magnetic  field  of  order  10 Tesla.  The  corresponding  power  produced  by  im3  of 
fluid  is  of  the  order  10 3 MW. 


6.2  The  Magnetohydrodynamics  of  Real  Generators 

The  magnetohydrodynamic  flow  in  real  generators  is  considerably  more  complex  than 
the  simple  flows  described  in  Parts  4  and  5.  MHD  generators  are  principally  configured  in 
linear  ducts  in  which  a  flow  proceeds  along  the  axis  of  the  duct  of  general  and  variable  cross 
sectional  shape  and  size.  The  Mach  numbers  of  interest  are  high  subsonic  or  moderately 
supersonic.  The  MHD  interaction  parameter  Su  will  be  of  the  order  6  — 10  for  commercial 
scale  (300MW(fc  —  2000AfW«fc)  systems.  The  typical  duct  aspect  ratio  will  be  of  the  order 
L/D  «  10.  The  viscous  Reynolds  number  will  be  of  the  order  J2e  =  106  for  typical  lm  scale 
cross  sections;  the  flow  will  therefore  be  turbulent  in  viscous  regions.  The  turbulent  MHD 
flow  in  such  ducts  will  exhibit  a  variety  of  complex  multi-dimensional  fluid  and  electrical 
phenomena  which  are  indicated  in  Fig.  6.3. 

6.2.1  Thermochemistry  and  Elei  rical  Conductivity 

The  single  most  important  parameter  for  a  combustion  driven  generator  is  the  elec¬ 
trical  conductivity  achievable  in  the  products  of  combustion  seeded  with  an  alkali  salt  of 
potassium  or  cesium.  This  quantity  (as  well  as  the  full  set  of  thermodynamic  and  transport 
properties  of  such  gas  mixtures)  is  well  predicted  by  detailed  thermochemical  calculations 
(Maxwell,  et.  al.,1973).  The  inputs  to  these  calcualtions  consist  of  the  heats  of  formation 
of  all  the  compounds  considered  in  the  mixture,  the  ionization  potentials  of  the  various 
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ions  (including  electron  afinities  for  negative  ions),  and  the  various  inter-particle  scattering 
cross  sections. 

The  Mollier  Diagram  generated  by  these  detailed  thermochemical  calculations  for 
typical  coal  combustion  products  is  shown  in  Fig  6.3  and  the  corresponding  electrical 
conductivity  is  shown  in  Fig  6.4.  It  can  be  seen  that  over  the  typical  expansion  envelope 
through  the  generator,  the  electrical  conductivity  ranges  from  105/m  to  1  S/m. 

6.2.2  Boundary  Layers 

The  fluid  in  the  anode  wall  layer  is  subjected  to  a  decelerating  axial  Lorentz  force 
which  is  comparable  to  that  in  the  core.  It  is  also  subjected  to  a  transverse  Lorentz  force 
which  is  set  up  by  the  Hall  current  induced  by  conductivity  non-uniformities  between  the 
boundary  layer  and  the  core.  This  force  (and  induced  Hall  current)  exist  even  in  a  Faraday 
connection.  The  forces  on  the  fluid  in  the  anode  boundary  layer  thus  tend  to  lift  the  fluid 
off  the  wall  and  to  induce  a  transverse  pressure  gradient.  This  transverse  body  force 
creates  an  anode  boundary  layer  which  grows  in  thickness  faster  than  a  conventional  layer, 
has  a  weaker  turbulence  production  within  it,  a  larger  boundary  layer  shape  factor,  and 
a  greater  tendency  to  stall.  Such  boundary  layers,  computed  with  the  multi-  dimensional 
equations  (2.3)  and  an  appropriate  turbulence  theory  (Demetriades,  Argyropoulos  and 
Lackner,  1971)  are  shown  in  Figs.  6.6-6. 8  for  a  generator  duct  flow  with  an  interaction 
parameter  Su  of  approximately  10. 

Like  the  anode  wall,  the  cathode  wall  layer  is  subjected  to  a  decelerating  Lorentz 
force  which  is  comparable  to  that  in  the  core.  The  nonuniformity  induced  transverse 
force,  however,  is  in  the  direction  to  drive  the  cathode  wall  layer  fluid  into  the  wall.  A 
corresponding  transverse  pressure  gradient  exists  as  shown  in  Fig.  6.6.  Contrary  to  the 
anode  wall  layer,  the  transverse  body  force  on  the  cathode  layer  intensifies  the  turbulence 
production  within  the  layer,  producing  a  thinner  boundary  layer  with  a  lower  shape  factor 
and  a  diminished  tendency  to  stall  compared  to  a  conventional  boundary  layer.  These 
general  features  of  the  cathode  boundary  layer  for  an  Su  **  10  scale  machine  may  be  seen 
in  Fig.  6.8. 

The  decelerating  Lorentz  force  which  the  fluid  in  the  sidewall  boundary  layer  expe¬ 
riences  is  significantly  diminished  compared  to  that  in  the  core.  In  a  uniform  conductiv¬ 
ity  flow,  the  velocity  profiles  in  the  boundary  regions  are  enhanced  by  the  Lorentz  force 
(Hartmann  effect)  as  described  in  Section  4.3.  When  conductivity  nonuniformities  are 
considered,  the  Hartmann  effect  becomes  more  pronounced  and  under  circumstances  of 
strong  interaction  (Su  >  1)  the  boundary  layers  possess  velocity  overshoots.  As  a  result, 
the  sidewall  boundary  layers  run  against  an  accelerating  pressure  gradient  compared  to  the 
core  leading  to  the  development  of  velocity  overshoots  and  negative  displacement  thick¬ 
nesses  (Fig.  6.7).  The  turbulence  production  in  the  part  of  the  boundary  layer  below  the 
maximum  velocity  point  is  greatly  intensified  over  that  in  the  anode  and  cathode  layers 
due  to  the  enhanced  shearing  rate  of  the  overshot  velocity  profile.  The  skin  friction  and 
heat  transfer  are  greatly  enhanced  on  the  sidewalls  (Fig.  6.8).  The  sidewall  layers  are  also 
free  of  Lorentz  forces  in  the  direct  normal  to  their  surfaces  since  this  is  the  magnetic  field 
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direction.  The  nature  of  the  sidewall  boundary  layers  in  a  300 MWth  machine  can  be  seen 
in  Fig.  6.7.  Boundary  layer  profiles  of  this  nature  have  been  measured  by  Daily,  Kruger, 
Self,  and  Eustis  (1976). 

In  a  conducting  sidewall  generator,  the  current  is  no  longer  constrained  to  enter  the 
electrode  across  the  magnetic  field  lines.  Under  some  circumstances  it  is  therefore  possible 
for  the  electrode  wall  boundary  layers  to  be  relieved  by  their  Lorentz  forces  as  the  current 
vector  aligns  itself  with  the  magnetic  field  in  the  boundary  layer  region  (J  x  B  =  0).  In 
such  a  situation  the  electrode  wall  layers  are  free  to  experience  the  analogous  acceleration 
which  the  sidewalls  experience  in  an  insulating  wall  generator  and  to  even  develop  velocity 
overshoots,  correspondingly  enhanced  turbulence  production,  shearing  rate,  and  enhanced 
heat  transfer  rates  (Markham,  Maxwell,  Demetriades.  and  Oliver,  1977.) 

The  actual  structure  of  the  flow  in  such  cases  will  be  complex  because  zones  of  en¬ 
hanced  Lorentz  force  will  develop  just  outside  the  boundary  layer  edge  as  the  current 
density  intensifies  but  still  contains  a  significant  component  in  the  transverse  direction. 
Further,  secondary  flows  which  transport  high  momentum  fluid  between  the  sidewalls, 
core,  and  electrode  walls  can  be  expected  to  modify  the  final  balance  of  momentum  in  the 
boundary  layer  regions.  High-interaction  conducting  sidewall  machines  therefore  require 
the  most  careful  analysis  of  multidimensional  and  secondary  flow  effects. 

6.2.3  Secondary  Flow  and  Magnetoaerothermal  Instability 

It  was  shown  in  Eq.  (2.11)  that  gradients  of  the  current  density  in  the  magnetic  field 
direction  are  a  source  of  vorticity.  When  such  current  gradients  arise  purely  from  velocity 
gradients  at  uniform  conductivity,  the  vorticity  production  is  negative,  i.e.,  the  vorticity  is 
damped  (Heywood.  1963).  On  the  other  hand,  when  the  current  gradient  is  induced  by  an 
electrical  conductivity  gradient,  the  vorticity  production  is  positive.  Such  nonuniformity 
induced  gradients  are  of  the  essence  of  MHD  induced  secondary  flow.  Because  of  the  low 
conductivity  region  near  the  electrodes,  a  Hall  current  Jx  is  induced  over  the  electrodes 
(Demetriades.  Argyropoulos  and  Casteel,  1970).  This  Hall  current  varies  in  the  magnetic 
field  direction  due  to  the  sidewall  cooling.  The  resulting  current  density  gradient  becomes 
a  source  of  axial  vorticity.  The  feedback  of  this  vorticity  on  the  conductivity  distribu¬ 
tion  can  lead  to  amplification  and  an  instability  known  as  the  magnetothermal  instability 
(Demetrides.  Oliver.  Swean  and  Maxwell,  1981). 

A  three-dimensional  evolving  flow  which  is  magnetoaerothermally  unstable  is  shown 
in  Figs  6.9  and  6.10.  The  case  considered  is  for  a  subsonic  segmented  Faraday  generator 
with  a  maximum  magnetic  field  strength  of  6 Tesla  and  an  interaction  parameter  Su  =  10. 
In  Fig.  6.9  the  streamlines  of  the  vortical  component  of  the  secondary  flow  velocity  field 
are  shown.  It  should  be  particularly  noted  that  even  after  only  lm  of  development,  MHD 
forces  have  already  dramatically  altered  the  conventional  turbulently  generated  corner 
secondary  flow  cells  transforming  them  to  fill  the  full  cross  section  of  the  channel.  In 
Fig.  6.10  the  inherently  magnetohydrodynamic  character  of  the  flow  is  evident  in  all  its 
characteristics  including  the  secondary  flow  effects.  In  particular,  the  growing,  intensifying 
current  concentration  in  the  central  anode  region  is  prominent.  This  magnetoaerothermal 
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current  concentration  has  brought  about  an  extreme  deceleration  of  the  primary  flow  in 
this  region  and  has  correspondingly  dramatically  reduced  the  local  skin  friction  in  the 
central  anode  region  and  increased  the  local  shape  factor  driving  this  region  of  the  anode 
boundary  layer  towards  separation  (Maxwell,  Swean,  Vetter.  Crouse,  Oliver,  Bangerter 
and  Demetriades.  1981). 

The  anode-sidewall  corner  cells  have  evolved  such  that  there  is  one  large  slowly  ro¬ 
tating  cell  which  tends  to  convect  hot  fluid  into  the  anode  and  relatively  colder  fluid  onto 
the  cathode  and  a  smaller  more  rapidly  rotating  cell  which  acts  across  the  anode.  In  the 
cathode-sidewall  corner  there  is  a  large  rapidly  rotating  cell  which  acts  along  the  sidewall 
boundary  layer  region  and  a  similar  but  smaller  cell  which  acts  over  part  of  the  cathode 
boundary  layer  region.  Together  these  cells  are  very  efficient  in  bringing  higher  tempera¬ 
ture  fluid  from  the  edge  of  the  boundary  layer  into  the  corner  region. 

Joule  heating  participates  with  secondary  flow  convection  in  determining  the  temper¬ 
ature  field.  The  locally  enhanced  temperature  and  conductivity  produce  regions  of  higher 
current  densities  which  in  turn  further  heat  the  fluid  through  Joule  heating.  This  is  the 
result  of  intensified  Joule  dissipation  in  certain  regions  which  is  not  convected  away  by 
the  secondary  flows  or  transferred  to  the  walls  by  turbulent  heat  conduction.  Typical 
temperature  distributions  in  such  a  situation  are  shown  in  Fig.  6.10b. 

6.2.4  Time  Dependent  Phenomena 

Time  dependent  phenomena  in  a  magnetohydrodynamic  flow  train.  Fig.  6.11.  may  be 
examined  with  the  quasi-  one-dimensional  system  of  Eqs.  (4.1)  (Oliver.  Crouse,  Maxwell 
and  Demetriades,  1980). A  basic  time  dependent  event  is  the  start-up  transient.  This 
dynamic  is  exhibited  in  the  space-time-  amplitude  plane  (Fig.  6.12)  including  the  ignition 
in  the  combustor  and  the  injection  of  the  ionizable  potassium  carbonate  which  takes  place 
approximately  8  milliseconds  after  ignition.  In  Fig.  6.13,  the  time  dependent  start-up  of 
the  flow  train  after  a  steady  oxidizer  cold  flow  has  been  established  is  shown. 

Of  particular  interest  are  the  extreme  excursions  of  the  electric  fields  during  this  start¬ 
up  transient.  The  axial  field  Ex  experiences  a  severe  overshoot  (Fig.  6.14)  because  the 
velocity  in  the  channel  is  high  before  the  seed  comes  on  and  the  Lorentz  force  decelerates 
the  flow  to  its  design  level. 


7.  MAGNETOHYDRODYNAMIC  ACCELERATORS 
7.1  Steady  Linear  Accelerators 

If  power  is  supplied  to  the  flow  and  Kv  >  1  as  described  in  Section  3.2.  the  Lorentz 
force  is  in  the  direction  of  flow  in  an  MHD  channel  and  the  working  fluid  is  accelerated. 
Steady  flow  channel  solutions  for  the  exit  state  of  the  fluid  in  terms  of  the  inlet  state  and 
the  applied  power  may  be  obtained  analogous  to  those  for  MHD  power  generators  (Jahn, 


Handbook  of  Fluid  Dynamics 


-  21  - 


©1987  STD  Research  Corporation 


1968).  The  velocity  enhancement  Ue/Uo  (where  Ut  is  the  exit  velocity  of  the  channel  and 
Uq  is  the  inlet  velocity)  is  a  function  of  the  interaction  parameter  Su  =  crB2L/pU: 


UJU0  = 


'  +  °<W)L 


(7.1) 


where  (  )  indicates  an  average  over  the  channel  and  the  pure  numbers  a,  a  have  differ¬ 
ent  order  unity  values  depending  upon  the  particular  assumptions  of  channel  operation 
(isothermal,  constant  field,  constant  current,  etc.)  A  pioneering  steady  flow  accelerator, 
the  Northrup  Mk  III,  achieved  a  thrust  level  of  approximately  AbNewtons  at  a  mass  flow 
of  1.3 g/s  of  argon  and  peak  current  densities  of  approximately  100 A/ cm2  (Demetriades  et. 
al.,  1963).  A  member  of  this  family  of  accelerators,  the  Mk  CA-1.  designed  for  accelerating 
air,  is  ilustrated  in  Fig  7.1. 


7.2  The  Magnetoplasmadynamic  Accelerator 

An  accelerator  of  considerable  interest  is  obtained  by  considering  a  steady  flow  ac¬ 
celerator  of  the  type  described  in  Section  7.1  in  cylindrical  coordinates  with  azimuthal 
symmetry  (Fig.  7.2).  For  this  configuration  of  anode  and  cathode  in  (r,0, z)  coordinates 
we  have  J  =  J(  Jr,0,0)  and  B  =  B(0,B^,0).  A  current  I  is  applied  between  the  anode  and 
the  cathode  and  the  Hall  current  is  suppressed.  The  magnetic  field  may  be  either  applied 
or  induced  by  the  self-current  ( Rm  >  1)  or  both.  With  only  variation  in  the  r  direction 
the  current  density  is  given  in  terms  of  the  applied  current  by 


Mr)  = 


2  xrL' 

The  induced  field  B«  is  determined  from  Eq.  (2.14b): 

dB9  _ 
dz 

and  thus  the  induced  magentic  field  is  given  by 


(7.2) 


(7.3) 


Be  = 


Ho  I 
2nr 


(1  -  z/L) . 


The  Lorentz  force  from  Eqs.  (2.14c)  and  (2.14d)  is  given  by  f  =  f(0,0,/z): 


(7.4). 


fz  =  JrB9.  (7.5) 

The  total  electromagnetic  thrust  force  acting  over  the  volume  defined  by  0  <  z  <  L,  rc  < 
r  <ra  is 


fx2nrdrd8dz  = 


4?r 


ln(r«/rc). 


(7.6) 
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The  thrust  force  is  thus  quadratic  in  the  applied  current.  For  typical  dimensions 
rc  «  1cm,  ra  «  10cm,  and  I  «  IkAmp,  a  thrust  force  of  the  order  1  /ANewton  is  generated. 
In  general,  the  current  pattern  will  possess  a  component  in  the  z  direction,  particularly  if 
the  anode  is  placed  downstream  of  the  cathode.  The  Lorentz  force  component  due  to  the 
presence  of  Jz  is  fr  =  —  JzBe  and  represents  a  focusing  (or  inward  pinching)  of  the  jet. 
When  the  Hall  effect  is  present,  the  aximuthal  symmetry  requires  that  the  Hall  field  Ee 
be  shorted  and  a  Hall  current  will  flow  given  by  Je  «  0Jr.  Correspondingly,  components 
Br ,  Bz  of  the  magnetic  field  will  be  induced  which  will  modify  the  focusing  and  thrust  forces 
and  also  lead  to  a  torque  rfe  =  r(JzBr  —  JrBz)  which  will  impart  angular  momentum  to 
the  fluid. 

The  MPD  accelerator  achieves  a  conducting  fluid  by  the  self-heating  of  the  working 
fluid  by  the  discharge  itself.  This  machine  operates  most  stably  at  low  pressure;  and  in 
this  regime  plasma  kinetic  effects  such  as  ionization  and  recombination,  charge  exchange, 
and  strongly  non-Maxwellian  behavior  play  an  important  role. 
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FIGURE  TITLES 


Fig.  4.1  Viscous  magnetohydrodynamic  flow  velocity  profiles  as  a  function  of  Hartman 
Number  Ha.  Lorentz  force  and  current  reverse  at  the  point  at  which  Ux  =  ( U ). 

Fig.  5.1  Magnetic  field  distribution  through  plasma  slab  moving  through  an  imposed 
external  magnetic  field  Bo  as  a  function  of  the  magnetic  Reynolds  number  Rm  =  pooU a. 
Note  that  for  Rm  >  10  the  field  (and  currents)  are  well  excluded  from  most  of  the  plasma. 
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Fig.  5.2  Planar  flow  at  finite  magnetic  Reynolds  number  Rm  =  1.  The  induced  magnetic 
field  isolevels  are  shown  In  (a),  the  electrodes  are  at  open  circuit  ( Hol/B .  =  0)  and  in  (b) 
the  electrodes  are  connected  to  a  load  ( Hol/B ,  =  1). 


Fig.  5.3  Planar  flow  at  finite  magnetic  Reynolds  number  Rm  =  7-  Induced  magnetic  field 
isolevels  are  shown.  In  (a),  the  electrodes  are  at  open  circuit  (/ XqI/B .  =  0)  and  in  (b)  the 
electrodes  are  connected  to  a  load  ( Hol/B .  =  1). 

Fig.  6.1  Magnetohydrodynamic  channel  configurations. 

Fig.  6.1a  Faraday  connected  linear  MHD  generator. 

Fig.  6.1b  Hall  connected  linear  MHD  Generator.  Faraday  circuits  are  shorted  and  power 
is  extracted  from  the  ends. 


Fig.  6.1c  Diagonal  conducting  wall  configuration.  Diagonal  conducting  bars  are  aligned 
with  design  point  equipotentials  and  power  is  extracted  from  the  ends. 


Fig.  6.2  Comparison  of  the  output  characteristics  of  diagonal  and  Faraday  generators. 


Fig.  6.3  Mollier  diagram  for  a  typical  coal  combustion  plasma  with  pre-heated,  oxygen- 
enriched  air,  and  potassium  carbonate  seed. 


Fig  6.4  Isopleths  of  scalar  electrical  conductivity  as  a  function  of  pressure  and  temperature 
for  a  typical  coal  combustion  plasma  with  pre-heated,  oxygen-  enriched  air,  and  potassium 
carbonate  seed. 


Fig.  6.5  Representation  of  multi-dimensional  fluid  and  electrical  phenomena  in  linear 
MHD  generator  ducts.  General  flow  field  involves  swirling  secondary  flow  resulting  from 
combustor  swirl  and  transverse  Lorentz  forces.  Current  distributions  are  concentrated  at 
electrode  edges  due  to  Hall  effects.  Heat  transfer  and  skin  friction  are  strongly  non-uniform 
due  to  current  concentrations  and  arcs.  In  slagging  generators,  hydrodynamic  slag  films 
interact  with  the  plasma. 


Fig.  6.6  Effect  of  Hall  current  on  pressure  distribution  between  electrodes.  Local  Hall  cur¬ 
rents  generate  local  transverse  Lorentz  forces  which  induce  transverse  pressure  gradients. 
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A  net  Hall  current  leads  to  a  pressure  on  the  cathode  wail  which  is  greater  than  the  anode 
wall  pressure. 


Fig.  6.7  Axial  velocity  profile  development  in  an  MHD  generator  of  scale  Su  =  10.  Sidewall 
velocity  overshoots  occur  because  of  cool  boundary  layer  regions  in  which  the  current,  and 
hence  the  Lorentz  force,  vanish. 


Fig.  6.8  Skin-friction  distributions  on  electrode  and  sidewalls  for  a  tvpicai  high  interaction 
MHD  generator  act  with  the  plasma.  Velocity  overshoots  on  sidewalls  lead  to  enhanced 
shear  and  skin  friction. 


Fig.  6.9  Development  of  the  streamfunction  formed  from  the  vortical  component  of  the 
secondary  flows  in  the  cross-section  at  lm  intervals  in  a  7m  long  MHD  generator  duct. 
Station  locations  are  relative  to  the  first  loaded  electrode. 


Fig.  6.10  Full  profile  structure  of  axial  velovity  and  vortical  secondary  flow  component 
streamlines  in  a  magnetoaerothermally  unstable  flow.  Transverse  Lorentz  forces  in  the 
core  resulting  from  the  Hall  effect  and  cool  walls  leads  to  the  development  of  a  set  of 
vortical  secondary  flow  cells.  These  cells  convect  hot  fluid  into  the  center  of  the  anode 
which  generates  a  further  increase  in  the  transverse  Lorentz  force  and  the  piasma. 


Fig.  6.11  MHD  generator  flow  train  with  control  points  which  induce  transient  effects. 


Fig.  6.12  Amplitude-space-time  representation  of  dynamic  variables  in  transient  flow  train 
response.  Any  variable  <p  has  an  amplitude  at  any  point  in  space  x  at  time  i  which  describes 
the  strength  of  the  disturbance.  The  projection  of  these  amplitudes  on  the  x  —  t  plane 
defines  a  trajectory  whose  slope  is  the  propagation  speed  of  any  disturbance  in  the  variable 


Fig.  6.13  Response  during  the  start-up  transient  of  a  large  (Su  «  10)  MHD  generator  flow 
train  as  manifested  in  the  pressure.  Initial  pressure  pulse  from  combustion  heat  release 
raises  the  combustor  pressure  to  a  peak  of  nealy  6afm  at  about  2ms  after  ignition.  This 
pulse  propagates  through  the  combustor  and  is  both  reflected  from  the  nozzle  throat  and 
transmitted  downstream  through  the  cold  oxidizer  flow.  The  downstream  wave  collides 
with  the  cold  flow  shock  leading  to  further  reflected  and  transmitted  waves.  The  second 
combustion  initiated  wave  propagates  downstream  in  hot  fluid  overtaking  and  merging 
with  the  first.  Hot,  conducting  fluid  propagates  through  the  system  beginning  at  7 ms 
after  ignition.  Once  conducting  fluid  enters  the  generator,  the  Lorentz  force  acts  upon  the 
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flow  and  the  original  cold  flow  shock  (which  exists  since  Lorentz  forces  are  not  present)  is 
dispersed.  Fluid  speed  is  U ;  speed  of  sound  is  c. 


Fig.  6.14  Hall  field  Ex  response  during  the  start-up  transient  of  a  large  MHD  flow  train. 
Fields  in  excess  of  4 kV/m  exist  for  a  period  of  about  10  —  ‘20ms  over  a  large  portion  of  the 
generator;  such  fields  can  result  in  axial  breakdown  arcs. 


Fig.  7.1  Cutaway  illustration  of  MI<  CA-1  accelerator  and  arc  jet  shown  mounted  on  thrust 
stand.  (Momentum  change  of  stream  is  measured  by  observing  thrust  change  on  strain 
gages.) 


Fig.  7.2  Cylindrical  geometry  MPD  accelerator. 
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